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ABSTRACT 

A major uncertainty in the dynamical history of the local group of galaxies originates from the unknown 
transverse speed of the Andromeda galaxy (M31) relative to the Milky Way. We show that the recent 
VLB A measurement of the proper motion of Andromeda's satellite, M33, severely constrains the possible 
values of M31's proper motion. The condition that M33's stellar disk will not be tidally disrupted 
by either M31 or the Milky Way over the past 10 billion years, favors a proper motion amplitude of 
100 ± 20 km s-i for M31 with the quadrant of a negative velocity component along Right Ascension 
and a positive component along Declination strongly ruled-out. This inference can be tested by future 
astrometric measurements with SIM, GAIA, or the SKA. Our results imply that the dark halos of 
Andromeda and the Milky Way will pass through each other within the next 5-10 billion years. 



1. INTRODUCTION 

The local group of galaxies provides the nearest labora- 
tory for the dynamics and evolution of galaxies. Its two 
dominant galaxies, the Milky Way and Andromeda (M31) 
are separated by a distance of ~ 700 kpc and are moving 
towards each other with a radial velocity of —117 km s^^ 
(Binney & Tremaine 1987). The unknown transverse ve- 
locity of M31 relative to the Milky Way introduces an 
ambiguity into the dynamical history of the local group 
and the likely dark-matter mass of its components (Pee- 
bles 1995; Peebles et al. 2001). The system potentially 
provides a local example for the dynamics that leads more 
generally to the formation of elliptical galaxies out of the 
merger of two disks (Barnes 1998). 

The estimated ages of stars within the disks of the Milky 
Way and Andromeda imply that these disks formed 10 
billion years ago, and their small scale height suggests that 
they had not experienced any substantial merger activity 
since then (Wyse 2000; Gilmore, Wyse, & Norris 2002; 
Peloso, da Silva, & Arani-Prado 2005; Toth & Ostriker 
1992). 

Recent VLB A observations of two H2O masers on oppo- 
site sides of Andromeda's satellite, M33, led to a geometric 
determination of its distance 730 ±168 kpc and total veloc- 
ity of 190 ± 59 km s~^ relative to the Milky Way (Brun- 
thaler et al. 2005). Unfortunately, no such masers were 
found in Andromeda, despite extensive search efforts (e.g. 
Imai et al. 2001). M33 has a thin H I/stellar disk that 
extends out to a radius of ~ 10 kpc (Corbelli & Schneider 
1997). Here we use the existence of this disk and its mea- 
sured position and bulk velocity in three-dimensions, to 
constrain the transverse velocity of Andromeda. In partic- 
ular, we use the constraints that: (i) M33's stellar disk was 
not disrupted by tidal interactions with the Milky Way and 
Andromeda over the past 10 billion years; and (ii) these 



galaxies were closer to each other more than 10 billion 
years ago as participants in the expansion of the universe. 
In §2 we describe our prescription for the mass profiles of 
the galaxy halos under consideration and in §3 we present 
our numerical results. Finally, §4 summarizes our main 
conclusions. 

2. HALO MASS PROFILES 

We would like to constrain the velocity vector of M31 
based on the condition that the stellar disk of M33 (as rep- 
resented by test particles on circular orbits at radii < 10 
kpc from the center of M33) was not disrupted during the 
past ~ 10 billion years. For simplicity, we assume a spher- 
ical distribution of enclosed mass, M(r), within each of 
the galaxy halos under consideration (M31, M33, and the 
Milky Way) and follow orbits of test particles in them. In 
order to solve for the orbits of these test particles we only 
need to specify the gravitational acceleration as a func- 
tion of radial distance r from the center of these halos, 
g = gr{r)er, where 

GM{r) v\t) 

9r = — :;3 — = — (1) 

with v{r) being the circular velocity as function of radius 
within these halos. 

Our test-particle approach ignores the effect of dynami- 
cal friction on the orbit of M33 inside the halo of M31. The 
additional effect of dynamical friction can only strengthen 
our constraints since it would bring M33 closer to the core 
of M31, where it would be more vulnerable to tidal dis- 
ruption. Hence, our simple tcst-particlc approach provides 
conservative constraints on the allowed range of possible 
orbits for the M33/M31 system. For an elaborate discus- 
sion on the impact of dynamical friction on the survival 
of satellite halos, see Taffoni et al. (2003) and references 
therein. 
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The simplest model of galaxy halos assumes a flat rota- 
tion curve, V = const, out to a truncation radius rj. Since 
all three galaxies under consideration have H I disks, the 
values of v can be inferred from their corresponding kine- 
matic data, namely ^ 110, 220, and 250 km s^^ for M33, 
the Milky Way, and M31, respectively. The truncation ra- 
dius is related to the total mass of the halo, M^, through 
the expression = GM^/v"^. For M33, rt is dictated by 
the tidal effect of M31. Given an orbit of M33, we find 
ri for M33 by searching for the largest radius at which a 
test particle initially on a circular orbit around M33 would 
remain bound to it. 

For M31 and the Milky Way, it is natural to associate 
rt with the virial radius, r^, inside of which the mean halo 
density is ^ 200 times larger than the mean density of the 
universe when these galaxies formed. More precisely, in 
a flat universe with density parameters in matter and 
f^A = 1 — Vim in a cosmological constant, the density of the 
halo in units of the critical density at the collapse redshift 
z of the halo is (see §2.3 in Barkana & Loeb 2001 for more 
details), 

Ac = IStt^ 82(i - 39(i^ (2) 
where d = — 1 and 

A halo of virial mass My that formed at a redshift z has 
a physical virial radius of 

= 84.5 X Miif[/(z)] -1/3 (^i±i^ ' /^-i kpc, (4) 

and a corresponding circular velocity at the virial radius, 

vy = 225.5 X MU^ifiz)]'/^ {H:^) ' • 

Here M12 = {M^jX^'^'^hr'^MQ), h is the present-day value 
of the Hubble constant in units of 100 km s~-'^Mpc~-'^, and 
f{z) = [(rim/r2^)(Ac/187r^)]. Our choice for the viriaUza- 
tion redshift of M31, M33, and the Milky Way is ^ ~ 1 
(Wyse 2000; Gilmore, Wyse, & Norris 2002). Nuclear dat- 
ing of the thin Galactic disk implies an age of ~ 8.3 ± 1.8 
Gyr corresponding to a formation redshift z ^ I (Peloso 
et al. 2005), as expected from other considerations for spi- 
ral galaxies (Hammer ct al. 2005). The total halo masses 
of the Milky Way and M31 were also constrained directly 
in the literature (Binney & Tremaine 1987; Gottesman, 
Hunter, & Boonyasait 2002; BattagUa ct al. 2005). We 
parametrize the masses of M31 and the Milky Way in 
units of 3.4 and 2.3 x 10^^ Mq respectively, and adopt the 
WMAP values for the cosmological parameters, namely 
Vim = 0.27, VIa = 1 - Vim = 0.73, and h = 0.71 (Spergel 
et al. 2003). 

To refine our constraints, we adopt a halo profile that 
follows from cosmological simulations of galaxy formation. 
Navarro, Frenk, & White (1997; hereafter NFW) found an 
analytic fit to results from their cosmological N-body sim- 
ulations of galaxy halos in a cold dark matter cosmology. 
The density profile in their model has an asymptotic ra- 
dial dependence of at large distances and at small 
distances. The model has two free parameters: one is the 
transition radius that separates these different power-law 

^ http://bima.astro.umd.edu/nemo/; see also http://www.ast.cam.ac. 



dependencies, r^, and the second is the so-called "concen- 
tration parameter" , c, which denotes the ratio between the 
virial radius and rg. 

The radial gravitational acceleration in the NFW model 
is given by 

GMyln{l + cx)-{cx)/{l + cx) 
r2 [ln(l + c)-c/(l + c)] ^' 

where x = r/vv with Vy given by equation (4). The 
present-day concentration parameter is correlated with the 
halo collapse redshift (Wechsler et al. 2002). 

3. NUMERICAL RESULTS 

Our goal is to examine a wide range of trial proper mo- 
tions for M31, compute the history of the M33/M31 /Milky 
Way system, and determine which proper motions lead to 
the survival of M33's stellar disk. Our numerical simula- 
tions used S. Aarseth's NBODYO code^, modified to allow 
the extended dark-matter potentials associated with the 
galaxies. M31 and the Milky Way were treated as ex- 
tended rigid bodies while the disk and halo of M33 were 
modelled as a collection of test particles moving in the 
combined, time-dependent potential of M31, the Milky 
Way, and M33. We used a Galactic Cartesian coordinate 
system in which the Milky Way is in the (x, y) plane, with 
the Galactic Center at (8,0,0) kpc, y is toward Galactic 
rotation, and z is toward the North Galactic Pole. 

The calculations were initialized to match today's posi- 
tions and motions of M33, M31, and the Milky Way. While 
we have a measured proper motion for M33, none exist for 
M31. Therefore, we assigned trial values for the proper 
motion for M31 in the eastward (a cos 5) and northward 
{5) directions and then converted them to the Galactic 
Cartesian frame. The observed radial velocity of M31 rel- 
ative to the center of the Milky Way is —117 km s"^, 
which converted to the Galactic Cartesian frame yields 
(56, —93, 43) km s^^. The radial and trial proper motions 
were combined to give the full velocity vector of M31 at 
the present epoch. For the distances of M31 and M33 from 
the Milky Way we adopt the respective values of 786 kpc 
and 794 kpc (McConnachie et al. 2004; 2005). 

M31 and the Milky Way were assigned total masses of 
3.4m/ and 2.3m/ x 10^^ M0 and truncation radii rt of 
235m/ and 207m f kpc respectively (consistently with their 
measured rotation curves for the isothermal sphere model). 
The scaling parameter m / was used to test the sensitivity 
of our results to the total dark matter masses (for recent 
empirical constraints, see Battaglia et al. (2005) and refer- 
ences therein). The dark matter distributions of M31 and 
the Milky Way were assigned an NFW halo mass profile 
with a concentration parameter c = 11.9 (Navarro, Frenk, 
& White 1997). Our results are not sensitive to reason- 
able variations in c. We also tested the truncated singular 
isothermal sphere as an alternative dark matter mass dis- 
tribution and found that it produced similar, albeit slightly 
greater, tidal stripping and heating of M33. 

The stellar disk of M33 at 10 or 5 Gyr ago was sim- 
ulated with typically 100 collisionless particles uniformly 
placed between 2.5 and 7.5 kpc of the center of M33 in 
the {x, y) plane. The outer cutoff simulates approximately 

k/sverre/web/pages/home.htm 
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Fig. 1. — Three examples of tidal stripping and heating of the stellar disk of M33. The locations of stars at the end of the N-body 
simulation are shown projected onto three orthogonal planes, as labeled at the top of the figure. These examples all used the full dark 
matter halos (mf = 1) and simulated histories over the past 10 Gyr. A 10 kpc scale is shown in the bottom left corner. Top: Simulation 
with (VacosSi Vs) = (100, —100) km s^-'-for M31, which resulted in essentially no tidal stripping (tidally stripped fraction of ~ 1%). Middle: 
Simulation with (VccosSi^s) = (—50,-50) km s~^for M31, which resulted in moderate tidal stripping (23%). Bottom: Simulation with 
(VacosSi^s) = (0)0) km s~^for M31, which resulted in heavy tidal stripping and heating (77%). 



the maximum extent of the stellar disk in M33 (CorbcUi 
& Schneider 1997). The inner cutolf was adopted for com- 
putational speed; since tidal heating and stripping starts 
from the outside and proceeds inward, the fate of test stars 
near the center of M33 is of little interest for our purposes. 

Since strong tidal stripping and heating of M33's halo 
is also a possibility, the dark matter distribution of M33 
was monitored by placing roughly 100 "mass-tracing" par- 
ticles at 1 kpc spacings out to the truncation radius of the 
galaxy. The enclosed mass profile for M33 was updated 
every 10^ yr from the locations of the mass-tracing par- 
ticles, azimuthally averaged and binned as a function of 
radius from the center of M33. As these mass-tracing par- 
ticles interacted with M31 and the Milky Way, some were 
tidally stripped, steepening the fall-off of mass with ra- 
dius in M33. The halo stripping and heating of M33 did 
not affect significantly its orbit (because of its low mass) 
but had some effect on the tidal disruption of its stellar 
disk. thcTc^orc started the calculation by integrating 
the coordinates of the three galaxies backwards in time 
to obtain initial conditions either 10 or 5 Gyr ago. This 
was implemented by reversing the sign of today's velocity 
vector for each galaxy. With the entire system specified 
in the past, the calculation proceeded forward in time to 
the present epoch and included the stripping and heating 
of M33's dark halo and stellar disk. 

Once the local group arrived at the present time, we 
displayed the locations of the "disk" stars in M33, as illus- 
trated in Figure 1. In some trials, the disk stars remained 
in orbit about the center of M33 with little change in ra- 
dial distribution or the original orbital plane, indicating 
little tidal stripping and heating. In other trials, many 



stars were moved to mtich greater radii, the mean orbital 
plane was rotated and warped, and the distribution was 
more triaxial than planar, indicating strong tidal heating 
after a few passages near M31 or the Milky Way. 

In order to evaluate quantitatively the degree of tidal 
stripping and heating, we produced a simple estimate of 
the percentage of stars "stripped" from M33. Stars moved 
to radii greater than 8 kpc or out of the initial stellar plane 
by more than 20% of their radial distance, were considered 
stripped. We determined the stellar plane in an approx- 
imate manner, by fitting straight lines to the projected 
distributions in the {x, z) and [y, z) planes, rotating the 
system back to the {x, y) plane, and calculating each star's 
distance from this plane (^;— distance). In general, three 
Euler angles are needed to correct the rotation of a system. 
However, for small tidal rotations of the stellar plane, our 
procedure is adequate since we start with the stars in a 
plane aligned with the Cartesian axes. For large rotations 
(> 30°), we optimized this procedure by considering a 7x 7 
grid of rotation angles spaced by 7° and centered on the 
initial results; for each test we calculated the percentage 
of stripped stars and adopted the rotations that produced 
the minimum stripping. 

Figure 2 displays the percentage of stars in M33 that 
were tidally stripped over a wide range of trial proper mo- 
tions for M31. For the galaxy masses listed above and 
a lookback time of 10 Gyr, large areas of proper-motion 
space for M31 result in significant (> 50%) tidal strip- 
ping of M33's stellar disk, including most proper motions 
smaller than ~ 75 km s^^. We repeated our calculations 
for two lookback times (5 or 10 Gyr) and for two estimates 
of total galaxy masses (given by m/ of 1.0 and 0.75). Re- 



4 




-100 100 




r. . I 1 

-100 100 




tj I I L_i 1 1 I L_i 1 1 I I 1 1 C 

-100 100 
M31 ^aooss (km s"^) 

Fig. 2. — The fraction of tidally stripped stars as a function of the trial proper motion of M31. Trial proper motions are in Equatorial 
coordinates (Right Ascension, a, and Declination, <5). Contours delineate 20% (light grey) and 50% (dark grey) of the total inunbcr of stars 
stripped in the simulation. Regions inside the dashed ellipses correspond to trial proper motions that satisfy the condition that the separation 
between M31 and the Milky Way was smaller 10 Gyr ago than it is today (the "timing" argument). Top: Simulations with full dark matter 
halos {mf = 1) and over the past 10 Gyr. The excluded region of proper motion is directed North- West. Middle: Simulations with reduced 
dark matter halos {mf = 0.75) and over the past 10 Gyr. Bottom: Simulations with full dark matter halos {mf = 1) and over the past 5 Gyr. 



diicing cither the lookback time or the masses of the galax- 
ies reduces the stripping, although significant stripping 
and heating still occurs over large areas of proper-motion 
space for M31. Since the real M33 does not appear to have 
undergone significant tidal heating, the darkly shaded ar- 
eas in Figure 2 would seem to be ruled out as possible 



proper motions for M31. 

We tested the sensitivity of our conclusions to uncertain- 
ties in M33's proper motion, the relative distances of M31 
and M33, and the initial orientation of the stellar disk in 
M33. In one test, shifting the motion of M33 by 70 km s~^ 
(randomly in the three Cartesian coordinates) increased 
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Fig. 3. — The separation between M33 and M31 as a function of lookback time for many trial proper motions of M31 (indicated the 
bounding axes). The scale for the individual plots is marked on the box near the lower left corner. The region inside the dashed ellipse 
corresponds to trial proper motions that satisfy the condition that the separation between M31 and the Milky Way was smaller 10 Gyr ago 
than it is today. All examples assume the full dark matter halos {nif = 1). 



the tidal stripping and heating by about 20%. Increasing 
(decreasing) the distance of M31 by 34 kpc [~ times 
the Icr precision of the observations in McConnachie et 
al. (2004,2005)] also increased (decreased) the stripping 
by ^ 20%. The conclusions of the paper would change 
significantly only if the distance to M31 is decreased by 
> 50 kpc (or through a similar increase in M33's distance, 
since the tidal effect depends primarily on the separation 
between M31 and M33). Finally, rotating the original stel- 
lar disk of M33 by 90° (from the {x, y) to the (x, z) plane) 
reduced the stripping and heating by about 30%. Overall, 
we conclude that our results are only sensitive at the level 
of tens of percent to uncertainties in the initial configura- 
tion of our simulations. 

One other constraint on the proper motion of M31 re- 
sults from the "timing" argument (Kahn & Woltjer 1959; 
Peebles et al. 2001). Because the Universe has expanded 
considerably over the past ~ 10 Gyr, galaxies (or proto- 
galaxies) should have been closer together in the distant 
past than they are today. Any large proper motion for M31 
violates this argument. Since our calculations produce a 
time history of the separation of M31 and the Milky Way, 
we plotted this separation as a function of time and de- 
termined the approximate range of M31 proper motions 
for which these galaxies were closer together 10 Gyr ago 
than today. The regions inside the dashed ellipses in Fig- 
ure 2 indicate approximately the proper motions of M31 
for which this occurs. For values of m/ < 0.5, we have 
found that M31 and the Milky Way were never closer to- 
gether in the past than today, for any of our trial proper 
motions of M31. Since this violates the timing argument, 
we do not show values of my < 0.5 in Figure 2. 

For the fiducial case of to/ = 1 and a lookback time 
of 10 Gyr (top panel), the allowed proper- motion param- 



eter space delineates an annular arc ^ 40 km s ^ wide 
that starts near Vm31 ~ (—90,0), goes counter clockwise 
through (0,-100) and (90,0), to (50,80) km s-^along 
Right Ascension and Declination. Additional small "is- 
lands" of moderate stripping and heating exist with 
smaller likelihood. 

Figure 3 shows the separation between M33 and M31 
as a function of lookback time for the various trial proper 
motions of M31. Since we initiate all orbits by integrat- 
ing M33's orbit back in time without stripping and then 
integrate these orbits forward in time with stripping, we 
lose time-reversal symmetry and so there is some varia- 
tion in the final separation between M31 and M33 among 
the different cases. Throughout the entire range allowed 
by the timing argument (delineated by the dashed line), 
M33 orbits around M31 several times as expected from a 
gravitationally-bound satellite. The repeated encounters, 
while not disruptive in the allowed range of proper mo- 
tions, might have still stripped material from the outskirts 
of M33 and could potentially account for the streams of H I 
gas (Braun & Thilker 2004; Thilker et al. 2004) and stars 
(Guhathakurta et al. 2005) that are observed in between 
the two galaxies today. 

4. CONCLUSIONS 

Figure 2 presents our limits on the proper motion of 
M31 based on the constraint that the thin stellar disk of 
M33 (simulated as an ensemble of collisionless test parti- 
cles on circular orbits) was not tidally disrupted over the 
past 10^° years through its orbit within the local group. 
This constraint and the condition that M31 and the Milky 
Way galaxies were closer in the past, favors a proper mo- 
tion amplitude of 100 ±20 km s~^ for M31 with the quad- 
rant of a negative velocity component along Right Ascen- 
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sion and a positive component along Declination strongly 
rulcd-out. The orbital angular momentum associated with 
the allowed range of proper motions implies that the dark 
halos of Andromeda and the Milky Way will pass through 
each other within the next 5 10 billion years. The main 
uncertainty in our results involves the dark halo masses of 
the local group members (compare the top two panels in 
Fig. 2). 

Our inference relies on the validity of various simplifying 
assumptions. For example, we simulated the dynamics of 
M33 without taking account of the full dark-matter distri- 
bution interior to the local group (including the mass in 
between the Milky Way and Andromeda). In addition, we 
have assumed that the tidal influence of external objects 
on the local group can be neglected, and we have ignored 



the addition of mass to the local group through accretion 
and infall of external matter. Obviously, the Milky Way 
and the Andromeda galaxies grew hierarchically with cos- 
mic time but the ages and low velocity dispersion of stars 
within their disks argue against major merging encounters 
over the past 10 billion years (Wyse 2000; Gilmore, Wyse, 
& Norris 2002). Our quantitative predictions provide an 
added incentive for future astrometric observatories such 
as SIM2, GAIA ^, or the SKA^, that will be able to test 
them (Shaya et al. 2003). 

We thank Andy Gould and Robert Braun for useful 
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